New Tools for Probing the Phase Space Structure 

of Dark Matter Halos 



Monica Valluri*, Victor P. Debattista f , Thomas Quinn** and Ben Moore* 

* University of Michigan 
^ University of Central Lancashire 
** University of Washington 
^University of Zurich 

Abstract. We summarize recent developments in the use of spectral methods for analyzing large 
numbers of orbits in A^-body simulations to obtain insights into the global phase space structure of 
dark matter halos. The fundamental frequencies of oscillation of orbits can be used to understand 
the physical mechanism by which the shapes of dark matter halos evolve in response to the growth 
of central baryonic components. Halos change shape primarily because individual orbits change 
their shapes adiabatically in response to the growth of a baryonic component, with those at small 
radii become preferentially rounder. Chaotic scattering of orbits occurs only when the central point 
mass is very compact and is equally effective for centrophobic long-axis tube orbits as it is for 
centrophilic box orbits. 
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INTRODUCTION 

While collisionless A/-body simulations produce dark matter halos that are triaxial or 
prolate, simulations which include the effects of baryons result in more spherical or 
axisymmetric halos[l, 2, 3]. It has been suggested that the change in shape could be 
the consequence of chaotic scattering of the box orbits that form the "back bone" of 
triaxial galaxies [4] . Alternatively, halo shapes might change due the adiabatic response 
of orbits to the change in the central potential. In a recent paper Valluri et al. [5] showed 
that by applying a spectral method to analyze large numbers of randomly selected orbits 
in A/-body simulations it was possible to clearly distinguish between these two options. 
In this paper we present a brief summary of some of their main results. 



FREQUENCY ANALYSIS OF N-BODY ORBITS 

Prolate and triaxial dark matter halos were formed from the merger of spherical NFW 
halos [6] and baryonic components (representing a disk, an elliptical galaxy or a massive 
compact nucleus) were grown adiabatically with time. The evolution of the halo was 
followed using PKDGRAV an efficient, multi- stepping, parallel tree code [7]. After the 
baryonic component was grown to full strength it was artifically "evaporated" to allows 
us to test for the importance of chaos in the evolution of the system [3, hereafter D08 



New Tools for Probing the Phase Space Structure of Dark Matter Halos 



February 4, 2010 



1 




"a *a 

FIGURE 1. For ^5000 pailicles frequencies before (co a ) versus after (C0b) the growth of a baryonic 
component are shown. Left: effect of an extended baryonic component (disk with 3 kpc scale length). 
Right: effect of a hard spherical point masses with 0. 1 kpc softening. 



and references therein]. In each system several thousand orbits were selected and their 
trajectories evolved in a frozen potential in each phase of the evolution of each halo. 
The initial triaxial/prolate phase is referred to as phase a, once the baryonic component 
is grown to full strength the halo is in phase b, and after the baryonic component had 
"evaporated" and the system returned to equilibrium the halo is in phase c. Each orbit 
was analyzed using a code that decompose its phase space trajectory to obtain its three 
fundamental frequencies of oscillation [8, 9]. The fundamental frequencies were then 
used to obtain a complete picture of the properties of individual orbits: to distinguish 
between regular and chaotic orbits, to classify regular orbits into major orbit families, to 
quantify the average shape of an orbit and relate its shape to the shape of the halo, and 
to identify the major resonant families of the halo which determine its structure [5]. 



RESULTS 



The largest of the three fundamental frequencies in each of the three phases of the 
evolution of our models (referred to as co a , (£>b an d co c respectively), can be compared to 
distinguish between adiabatic and chaotic evolution. In the case of a primarily adiabatic 
response, particles deep in the potential (large (O a ) are expected to experience a greater 
increase in frequency than particles further from the center. In Figure 1 we compare 
the frequencies of orbits in an initially triaxial halo (co a ) with their frequency in the 
presence of a baryonic component ((Ob). For the extended disk (Left) (as well as other 
extended baryonic components), (Ob increases monotonically with (O a with fairly small 
scatter, indicating that the orbits in this potential responded adiabatically (the dashed line 
shows the 1:1 correlation between the two frequencies). However, a hard central point 
mass (Right) results in significant scattering in a>& with small values of (O a (i.e. weakly 
bound orbits) having some of the largest values of CO/,. (Ob sometimes decreased instead 
of increasing providing further evidence for chaotic scattering in this case. 

The orbital frequencies were used to distinguish between regular and chaotic orbits 
and to classify the regular orbits into major families. We studied both triaxial and prolate 
halos and our orbital analysis showed that the initial triaxial halo was composed of 84- 
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FIGURE 2. The change in frequency of an orbit (0 ac versus pericentric radii r pe ri for triaxial halo (Left) 
and prolate halo (Right). 



86% box orbits, 11-12% long-axis tube orbits, 2% short axis tubes, and 1-2% chaotic 
orbits. In contrast the prolate halo had 15% box orbits, 78% long-axis tubes, 7% short- 
axis tubes and no chaotic orbits. To test the hypothesis that it is the box orbits that are 
most significantly scattered by a central point mass we define the fractional change Aco ac 
in the frequency of an orbit from its value in the initial halo (co a ) to its value after the 
baryonic component was evaporated {(O c ) to measure the amount of chaotic scattering. 

Figure 2 shows that orbits with smaller pericenter radii r per i experience a significantly 
larger change in frequency A(0 ac than orbits at large pericenter radii. This is true for both 
the triaxial halo (Left ) which is dominated by centrophilic box orbits and the prolate 
halo (Right) which is dominated by long-axis tubes. Thus chaotic scattering is equally 
strong for the centrophobic long-axis tube orbits and centrophilic box orbits contrary to 
the prevailing view. 

In a triaxial potential in which the major, intermediate and short axes are along the 
Cartesian directions x,y,z respectively, the oscillation frequencies satisfy the condition 
| CO* | < \(0y\ < \(O z \ for any orbit with the same over-all shape as the density distribution. 
For such an orbit we can define a shape parameter % s = \(O y /(O z \ — \(O x /(O z \ which is 
positive for orbits with elongation along the major axis of the figure, with larger values 
of Xs implying a greater degree of elongation. Since orbits closer to the central potential 
are more significantly affected by the baryonic component, we expect them to become 
rounder (% s — > 0) than orbits further out. Figure 3 shows how the orbital shape parameter 
Xs for various orbital types varies as a function of the pericenter radius r per i for a triaxial 
halo (Left) and after a disk was grown in this halo with symmetry axis parallel to the 
short axis (Right). In each plot the curves show the mean shape distribution of orbits of 
a given orbital type at each radius, with orbital types indicated in the line-legends. In the 
initial triaxial halo boxes, long-axis tubes and chaotic orbits are significantly elongated 
(Xs ~ 0.35) both at small and large radii. After the growth of the short-axis disk the 
orbits at small radii become axisymmetric (x s ~ 0) while orbits at large radii (especially 
boxes) remain quite elongated. This change in orbital shape at small radii was seen in 
all halos regardless of the radial scale length of the baryonic component. 
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FIGURE 3. The mean orbital shape parameter % s for different orbital types (indicated by line legends) 
vs. pericentric radius r pel j. 



IMPLICATIONS 

The analysis of fundamental frequencies of orbits in Af-body halos is a powerful tech- 
nique that allows us to identify the primary physical processes that cause halo shapes to 
change in response to the growth of a baryonic component. We confirmed the conclusion 
reached by D08 that chaos is not an important driver of shape evolution but found that 
significant chaotic scattering does occur when the baryonic component is in the form of 
a hard central point mass (of scale length ~ 0.1 kpc). Regardless of the original orbital 
composition of the triaxial or prolate halo, and regardless of the shape and radial scale 
length of the baryonic component, halos become more oblate because orbits closer to 
the center of the potential become axisymmetric. Although the resultant halos are al- 
most oblate, their orbit populations contain orbits (boxes and long-axis tubes) that are 
characteristic of a triaxial rather than axisymmetric potential. 
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